基于有限体积虚拟区域方法的两相流动直接模拟

为提高计算效率，提出有限体积法离散下的虚拟区域颗粒两相流动直接模拟方法.在控制方程中加入相应的虚拟区域源项，保证了颗粒内部的刚体运动特性.该源项中含有颗粒信息部分及流体信息部分.在每次迭代后，对源项中的流体信息部分进行更新，从而更好地保证颗粒内速度的刚体分布.计算静止颗粒圆柱绕流及单个颗粒的沉降过程，验证了算法的准确性.     

A new formulation for fictitious mass of viscous dynamic relaxation method

Abstract

In this article, a new relationship is proposed for the fictitious mass of viscous dynamic relaxation (DR) method. First, incremental equations are derived for DR steps. Using transformed Gershgörin theory, a new relationship is achieved for fictitious mass of viscous DR by formulating modified time step ratio. This procedure presents a new algorithm for the viscous DR method. To evaluate the numerical efficiency of the proposed method, some 2D and 3D truss and frame structures are analyzed with elastic linear and geometrically nonlinear behaviors. Results show that by using the proposed algorithm for fictitious mass, the convergence rate of the viscous DR method is improved so that the proposed algorithm presents the structural response with lower iterations in comparison with other common DR techniques.

Communicated by Joerg Fehr.     

Infrared thermographic evaluation of large composite grid parts subjected to axial loading

Many composite parts, such as laminated panels and grid-like shells, operate under high mechanical loading. Evaluation of their structural integrity is crucial to ensure the long-lasting operation of critical components. Since testing a structure under full or “proof” load might be dangerous for personnel, it would be preferable to use a remote, rapid inspection technique. This paper describes a practical application of IR thermography to the inspection of large composite parts used in the aerospace industry. This work has used just one cycle of increasing load from zero load to failure, and this was done for both for tensile and compressive loads. It is shown that, during the formation of micro-defects in polymeric composites, about 40 % of the total dissipated energy is expended for material heating, while about 60 % is related to material damage accompanied by an increase in the defect concentration. Non-uniform composite deformation causes temperature anomalies, whose amplitude may reach 1.5–2.5 °C at a load of about 50–60 % of the limit load.     

Temperature-dependent ratcheting of PTFE gaskets under cyclic compressive loads with small stress amplitude

Uniaxial stress-controlled ratcheting experiments on PTFE gaskets under cyclic compressive loads with small stress amplitude were performed. The effect of temperature on the deformation behavior was considered. Results showed that the compressive modulus decreases rapidly when the temperature increases from 100 °C to 200 °C. Compressive ratcheting deformation with cycles increase significantly with the increases of temperature. The ratcheting deformations at 100 °C, 150 °C and 200 °C are nearly two, three and five times that at room temperature, respectively. Most of ratcheting deformation mainly occurs during the first 20 cycles because the subsequent ratcheting rate and strain range are small and much less than those in the previous cycles. The accumulated deformation under cyclic loads with small stress amplitude is relatively approach to the static compressive creep with the same peak stress. Therefore, the accumulated deformation with time of PTFE gaskets obtained by cyclic compression with small stress amplitude can be estimated by the corresponding static creep deformation with good accuracy under the approximate stress rate and the same temperature, especially at room temperature.     

Energy release rates in a constrained epoxy disc with Hookean and Mooney–Rivlin materials

Assuming that the disc material can be modeled either as Mooney–Rivlin or as Hookean and the steel ring enclosing the disc as Hookean, the energy release rates as a function of the crack length are evaluated and compared. Two loadings are considered––one in which the surface of the star shape hole in the disc is loaded by a uniform pressure and the other in which the temperature of the composite body is uniformly raised. It is found that the linear and the nonlinear analyses give qualitatively similar results for the two loadings. For each load, the energy release rate increases with an increase in the starter crack length, reaches a maximum value and then decreases gradually.     

On Determining the Percolation Reynolds Number and the Characteristic Linear Dimension for Ideal and Fictitious Porous Media

Various versions of representations of the percolation Reynolds number for porous media with isotropic and anisotropic flow properties are considered. The formulas are derived and the variants are analyzed with reference to model porous media with a periodic microstructure formed by systems of capillaries and packings consisting of spheres of constant diameter (ideal and fictitious porous media, respectively). A generalization of the Kozeny formula is given for determining the capillary diameter in an ideal porous medium equivalent to a fictitious medium with respect to permeability and porosity and it is shown that the capillary diameter is nonuniquely determined. Relations for recalculating values of the Reynolds number determined by means of formulas proposed earlier are given and it is shown that taking the microstructure of porous media into account, as proposed in [1, 2], makes it possible to explain the large scatter of the numerical values of the Reynolds number in processing the experimental data.     

Stability of longitudinally corrugated cylindrical shells under uniform surface pressure

The buckling problem for longitudinally corrugated cylindrical shells under external pressure is solved. The solution makes practically exact allowance for the geometry and buckling modes of the shell. The inaccuracy of the results is due to the assumption that the subcritical state is momentless. Shells consisting of cylindrical panels of smaller radius and noncircular shells with sinusoidal corrugations are analyzed for stability. The practical applicability of such shells is demonstrated __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 11, pp. 66–79, October 2007.     

Application of a kinematic hardening viscoplasticity model with thresholds to the residual stress relaxation

The life analysis of engine components needs to take into account the residual stress relaxation induced by cyclic service loads. The paper recalls a new class of constitutive equations for cyclic viscoplasticity, using a series of kinematic hardening models with thresholds. The equations are introduced within a recently enlarged thermodynamic framework. Some attention is focused to the relations with multisurface approaches and to a specific determination procedure of the model parameters. The new model is applied to the calculation of the near surface residual stress relaxation after shot peening, when the structure is submitted to cyclic service loads. The simulated stabilized residual stresses are in good accordance with experimental results obtained on an N18 disk alloy at 650°C. In comparison, the classical model without threshold predicts the complete vanishing of the residual stresses, which is not satisfactory.     

Short-term microdamage of a physically nonlinear particulate material under a combination of normal and tangential loads

The structural theory of short-term damage is generalized to the case where the undamaged components of a particulate composite deform nonlinearly under loads that induce a compound stress state. The basis for this generalization is the stochastic elasticity equations for a particulate composite with porous components whose skeletons deform nonlinearly. Damage in a microvolume of the material is assumed to occur in accordance with the Huber-Mises failure criterion. Balance equations for damaged microvolume are derived for the physically nonlinear materials of the components. Together with the macrostress-macrostrain relationship for a particulate composite with porous nonlinear components, they constitute a closed-form system of equations. This system describes the coupled processes of physically nonlinear deformation and microdamage. Algorithms for calculating the microdamage-macrostrain relationship and plotting stress-strain curves are proposed. Such curves are plotted for the case where the composite is subjected to a combination of normal and tangential loads, and microdamages occur in the linearly hardened matrix and do not in the linearly elastic inclusions. The stress-strain curves are examined depending on the volume fraction of inclusions and presence of tangential stresses __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 12, pp. 48–57, December, 2006.     

Flexural-torsional bifurcations of a cantilever beam under potential and circulatory forces I: Non-linear model and stability analysis

The stability of a cantilever elastic beam with rectangular cross-section under the action of a follower tangential force and a bending conservative couple at the free end is analyzed. The beam is herein modeled as a non-linear Cosserat rod model. Non-linear, partial integro-differential equations of motion are derived expanded up to cubic terms in the transversal displacement and torsional angle of the beam. The linear stability of the trivial equilibrium is studied, revealing the existence of buckling, flutter and double-zero critical points. Interaction between conservative and non-conservative loads with respect to the stability problem is discussed. The critical spectral properties are derived and the corresponding critical eigenspace is evaluated.